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The  Eleciroii  Nndeir  Dynamics  0^^)  snffoadi  is  developed  for  a  linear  chain  in  a 
parametrized  model  inqtiied  by  the  FFP  (Pariser-Farr-Pople)  model  Particular  attention  is  given 
to  die  model  parameters,  and  die  dicnoe  of  basis  funcdcms  in  diis  time-dqiendent  theory.  The 
resnltiiii  equmioos  of  motion  indude  dectromc-vibrational  coiqilings.  Erplidt  analysis  of  the 
sinqdest  nmdd  leads  to  cotqdmg  between  die  highest  frequency  longhudi^  vitnadonal  mode 
and  die  dectroos. 


^4“*  15943 

mrnm 


94  5  26  105 


r.  Boa  SM.  s-nuo  OnMlB.  S«9dM 


DTIC  QUALITY  USSPEOfBD  1 


report  documentation  page 

'Apw  Aflgfovte 

OA«a  Na.  OfOd-OtiO 

i7iiiiR¥uiraiYiu^ owMij  i.  Mpcrn  oAti  imroarvYH  AM 

May  16,  1994  Tachnieal 

9  OATIS  COVIAIO 

leport 

A.  Tim  ANo  sustim  « 

A  Modal  For  Bleetron  Muclaar  Coupling 

i 

S.  AUMOMIG  NUMStAS 

N00014-93- 1-0122 
BT4131072 

4.  AU^HOAtS)  j 

J.'L.  Calais,  B.  Dauaana,  and  Y.  Obvn  ^ 

1 

f.  naFOMNMQ  oaoANaAnoN  namks)  ano  aooacsmis) 

Unlvsralty  of  Florida 

Quantun  Theory  Project,  HM  362 

Gainesville,  FL  32611-8435 

A  ASATOAMING  OAGANUATION 

AfPOAT  NUMSIA 

Tech.  Rep.  No.  3 

t.  SPONSOAMa.  lllbM'fdAIN^  A<icNO  NAMKS)  ANO  AOOAISKiS) 

Office  of  Maval  Bosearch 

Chenlatry  Prograa 

800.  Mortb  Quincy  8t. 

Arlington,  VA  22217-5000 

10.  SPONSOAiNO/  M0NIT0AIN6 

AGiaCV  MPOAT  NUMSIA 

11.  SUAPUMSMTAAV  NOTIS 

Subaltted  to  J.  Chea.  Phys.  (1994) 

12a.  eaTfMUTtON/AVAAASajTV  STATtMSNT 

This  docuaent  haa  been  approved  for  public  release  and 
ealei  Its  distribution  la  unllalted 

12fe.  OiSTAlSUTION  coot 

Unliad.ted 

1A  AasnuCT  (mmmim  200  ¥tCfau 

The  Electron  Nedear  Dynamics  (END)  approach  is  devdoped  for  a  linear  chain  in  a 
paiamettiaed  model  inspiied  by  the  PFP  (Pariser-Pair-Pople)  model  Paiticnlar  attendon  is  given 
to  die  ffwIH  parameieiSt  and  die  dioice  oi  basis  limcdoos  in  this  time-dependent  dieoiy.  The 
lesnlting  eQioatioos  of  motion  inclode  electxonio-vibratiooal  couplings.  Explicit  analysis  of  the 
yitlipittt  model  lf*t«  to  ctmpling  between  the  highest  fieqoency  longitudinal  vitnational  mode 
and  the  dectraos. 
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Tbe  raoendy  developed  Etoctnm  Nudear  I^namics  (END)  method  [1-10]  ofifen  {wcnniaing 
opportooities  toweidt  an  integrated  undeiftanding  of  electrcmic  and  nuclear  motion  and  their 
ooopling.  Instead  ot  die  tiacBtumal  sqMratum  of  dial  problem  into  basically  three  stqts  —  the 
calculation  ci  a  potential  energy  surface,  the  dmennination  of  an  analytical  n^nesentatum  of  that 
surfitee,  and  die  dynamics  on  the  sutfiaoe  —  the  END  method  attacks  die  problmn  in  one  stq>. 
This  has  bodi  concqicoal  and  computatitmal  advantages.  The  general  fiam^oik  of  die  END  is 
a  vaiiadoiial  appcoadi  towards  die  solutitm  ci  the  time-dqiendent  SchrOdinger  equatimi  for  all 
die  partides  of  interest  in  die  system. 


Bn^ldoying  die  time-dependent  variational  princ4>le  (TDVP)  [11,  12],  a  set  of  first-order 
coiqded  d^eteatial  equations  for  parameters  duutacterizing  a  time-dependmit  state  vectm  is 
derived.  These  time  dqiendent  parameters  define  an  approximate  solutitm  of  the  time  dependent 
SchrOdinger  equation  for  the  system  under  study.  The  parameters  can  be  of  "electronic”  type, 
sndi  as  coefBcksits  of  basis  fimctioos  in  molecular  or  crystal  mbitals,  or  of  "nuclear”  type,  Le. 
diaracterizing  die  stale  vectors  for  the  nudd.  The  meduxl  has  die  flexibility  to  accommodate 
ddier  a  quantum  mechanical  or  a  classical  treatment  oi  die  nucleL  When  the  nudd  are  treated 
as  classical  particles  die  niidear  coordinaies  and  momenta  are  the  parameters  that  evolve  in  time. 

The  END  dieoiy  can  be  viewed  as  a  geoenl  iq^NPoadi  to  study  time  evolnthm,  and  thus  also, 
duon^  appropriate  Fourier  transforms,  qiectnt  in  general  scattering  or  bound  state  problems. 
A  particular  choke  of  die  form  of  stale  vector  or  wavefhnction  and  associated  basis  sets  yields 
a  certain  leaUzatioii  of  END  or  in  odier  wmds  a  particular  modd  within  die  END  framework. 
When  parameters  are  indqiendent  of  time  die  equations  of  motion  reduce  to  die  extmnum 
conditions  of  die  time-indqiendeot  variationd  prind^  The  computer  program  ENDyne  [13], 
tdiich  inqiiements  dte  dieory,  can  dins  also  be  vaod  for  geometry  optimization  calculations.  Widi 
frozen  nudd  the  dynamical  equations  can  accomplish  wavefunction  optimization  and  in  general 
^ultaneoQS  geometry  and  wavefhnctkm  optimization. 

The  simidest  possible  approximation  of  END  dieory  invdves  dasskal  nudd  and  a  single 
determinantal  waveAmction  for  die  ekctrons.  The  latter  choke  nu^  seem  scmiewhat  restricted, 
however,  widi  appropriate  dioioe  of  wavefunction  parameters  the  determiiumtal  wavefuncticm 
can  be  given  die  form  of  a  cohereitf  sane.  In  this  so  called  Thouless  form  [14]  of  die  dderminant 
the  parameters  that  determine  die  degree  of  mixing  of  the  basis  fimetions  are  diosen  in  sndi  a 
wqr  that  as  diey  diange  in  time  all  possible  determinantal  states  in  die  given  basis  can  in  princi|de 
be  acoesaed.  The  dynamical  qnn  orbitals  dun  make  iqi  die  Thouless  determinant  are  comj^ 
nonordiogooal  functions,  whidi  can  permit  general  qnn  orbitals  of  various  kinds  in  line  with 
die  work  of  Pulrattome  [15].  The  computer  code  ENDyne  is  implemented  at  this  simplest  levd 
of  approximarion  However,  the  dieory  has  been  worl^  out  for  a  general  multictmfigmational 
state  vector  for  die  dectroos  [6]  and  for  a  wave  packet  treatment  of  die  nuclei  [5]. 

Mote  of  die  ^i|dicatians  of  END  carried  ote  so  far  are  concerned  with  reactive  ccdliskn 
problems  invdving  small  qiecies,  and  in  particular  dectron  transfisr  reactions  for  such  systems. 
Modifications  necessary  for  larger  systems  induding  extended  pdymerk  systems  are  nndte  way. 
The  preaem  pi^  employs  the  END  to  develop  a  cemsistent  procedure  to  describe  the  coiqiliitg 
bttwoea  dectnoic  and  nuclear  motion  in  extended  systems.  This  will  bring  two  areas  together 
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wUdi  h«ve  $0  ftr  often  been  treated  aqwfately,  namely  lattior.  dynamics  and  band  theofy.  Only 
anaDi  vibcaions  oi  die  nuclei  (ion  cores)  aroond  dieir  equilibrium  poritkms  are  considered  in 
dm  Tnodri 

The  pmpose  of  dm  present  p^per  is  to  identify  in  some  detail  and  address  a  numba  of 
problems  cofmected  with  dm  iqiplicatioo  of  dm  END  mediod  to  an  extended  system.  The  long 
lange  aqmct  of  dm  Coulomb  interactioo  becomes  a  pardcnlariy  important  issue  in  dm  treatmem 
of  extended  systems,  with  wavefuncdons  diat  must  extend  duoughout  dm  system  [IQ.  Both  in 
traditional  laftioe  tfynamics  and  band  dmocy  dm  associatfd  problems  are  w^-4nown  [17],  and 
it  is  essential  toke^  dmm  In  mind  srimn  develo|dng  anew  approach. 

Instead  of  attempting  directly  a  full  nb  biido  treatment  a  model  system  is  discussed  vrimre 
some  of  dm  siqm  in  dm  development  can  be  handled  analytically.  Urn  Pariser-Pair-Eople  (FFP) 
model  [18,  19]  was  originally  developed  for  incorpocatiog  dectron-dectron  interaction  in  dm 
treatment  of  x-electron  mdlecolar  systems.  Dm  Hamiltonian  [20],  vriiicb  has  grown  out 
of  dm  original  wodc,  has  been  iq^Ued  to  a  wide  range  of  many-electron  problems.  It  has  beat 
qqdied  at  several  levels  of  approodmadon  not  dm  least  to  dm  treatment  d  pdymers  (see  for 
instanoe  [21]  and  [22]  and  references  therein).  This  PPP  Hamiltonian  forms  dm  starting  point  of 
dm  present  p^mr.  The  particular  features  d  dm  END  mediod  and  dm  fact  diat  it  coadders  bodi 
nuclear  and  dectronic  dynamics  requires  a  careful  analysis  and  qmcification  of  the  traditional 
parameters.  An  essential  aqma  tte  PPP  modd  is  that  the  parameters  can  be  chosen  so  as  to 
preserve  dm  long  range  properties  of  dm  Coulomb  interaction.  The  energy  ttf  dm  system  for  a 
stttic  lattice  in  diis  modd  vrill  therefore  be  propottiond  to  dm  ruimber  of  atoms  (dm  volume)  of 
dm  qrstem  and  it  is  important  to  study  how  dm  energy  depends  on  dm  small  nuclear  vibrations 
around  dm  equilibtinm  poritions.  A  harmonic  mproxlmtfion  in  terms  of  dmse  deviations  is 
dmrefore  natural  A  harmonic  q^roximatton  of  the  told  energy  is  adopted  widi  reqmct  to  dm 
electronic  parameters.  The  modd  dmn  rqxesents  a  idly  coupled  dectronic>noclear  problnn 
at  dm  levd  of  dm  Random  Phase  Approximation  (RPA)  or  Immuiaed  DmO’Dqmndent  Hartree 
Pbck  (DHIF).  This  qiproadi  has  beia  discussed  in  generd  terms  in  reference  [3]. 

Dm  main  part  of  dm  paper  is  planned  as  follows.  Dm  basic  END  ftameworit  is  summarized 
indmnextaectikm.  hi  section  IE  a  number  of  probtons  associated  vrith  the  treatment  of  extended 
^sterns  ate  reviewed,  widi  particular  attemion  given  to  dm  aqwration  of  long  and  shmt  range 
forces  and  also  to  dm  importance  of  working  widi  quantities  dun  scale  coriecdy  with  dm  size 
of  dm  system,  hi  section  IV  the  particulars  the  chosen  modd  are  discussed.  Dm  question  of 
how  to  diooae  dm  electiooic  basis  fbnctimis  for  dm  END  treatment  is  discussed  in  a  sqmrate 
aectioiL  hi  sections  VI,  VE,  and  VIE  explidt  appUcations  of  END  to  duee  cases  are  carried  out 
Rrst,  in  order  to  commct  to  more  traditional  treatments,  dm  cases  with  mily  nuclear  parameters 
(section  VI),  and  only  dectronic  parameters  (section  VE)  are  treated.  Dmn  in  section  VIE  the 
general  case  is  ooDStdered  with  bodi  electronic  and  nuclear  parameters.  RnaOy,  the  results  are 
summarized  and  further  possible  extensions  are  discussed  in  the  last  aectitm. 

U.  Thft  END  equations  for  a  aingie  determinant  and  classical  nuclei. 

A  lattice,  periodic  in  one  dimensicm,  is  ctmsidered,  that  cmisists  of  M  unit  cells  of  length 
a  vridi  one  atom  per  unit  oelL  <3eneralization  to  an  arbitrary  number  of  aUnns  per  unit  odl 
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it  itnigiitforwtid.  Periodic  boundaiy  conditioot  tie  imposed  so  dist  ill  wivefiincti(xis  satisfy 

d  « 

=  (HI) 

Hie  im^erm  labels  the  H/ cells  in  die  Boni*voQ  Kinnin  region  (BK)  ~Af/2  <  m  <  Af/2  - 1. 

Ihe  actual  nuclear  (core)  positicns  are  denoted 

^m)  SB  111  +  ^m)  (IL2) 

udwie  ril  is  a  point  in  cdl  m  defined  by  the  equililxium  positicm  for  that  nucleus.  In  other  woids 
tiiis  is  die  anchor ptrifit  of  die  nudeus  (core)  and  ^m)  is  die  displacement  from  the  equilibrium 
of  diat  nuclens  (core).  The  presentation  is  limited  to  small  vibrations  around  the  equilibrium 
positions  so  diat  p  «  R.  The  nucleus  in  cell  m  has  momentum  P{m).  The  diqilaoements 
and  therefore  the  nuclear  positions,  as  wdl  as  die  momenta  depend  on  die  time  parameter 
t,  twUle  die  andior  points  rh  are  time  indqiendent 

The  ekctrons  in  the  system  are  described  1^  a  single  determinaiit  built  from  dynamic  spin 
orbitals  ss  (f;^)  with  (  the  qnn  cootdinree) 

K 

w«)  =  ^i«)+  E  (03) 

i-jr+1 

The  rank  ci  die  qpin  orbital  basis  is  K  and  the  first  AT  of  the  basis  spin  cnbitals  make  iq> 
die  so  called  refptence  determinant  The  time  dependent  electnmic  wavefunction  is  defined  as 

l*>  *  <iet{x<(ei)}  (n.4) 

rriiere  a  is  used  as  a  collective  symbd  for  all  the  N(K~N)  complex  paiametea  z^i  s  zji(t).  This 
form  of  die  determinantal  wavefimctimi  has  die  form  of  a  generalised  coherent  state  [23, 24, 3]. 
It  is  mmoimalized  and  die  normalisation  is 

5(**, »)  s=  (x|s)  «  det{l  +  »*»}.  (IL5) 

Oassical  treatment  of  the  nuclei  can  be  ctmsidered  as  die  narrow  wavepacket  limit  of  a 
product  of  frozen  Oanssians,  The  TDVP  is  used  to  arrive  at  a  s^  of  coqiled  differential  equations 
for  die  efectron  nuclear  dynamics  by  making  stationary  the  quantum  mechanical  actimi  [8].  The 
dynamical  equations  can  be  mqnereed  as 

(iC  0  0  0\/i\  /dBldar\ 

0  -iC*  0  0  I  i*  1 

0  0  0  -1  I  R  “  I  dE/dR  I 

0  0  1  0/\P/  \aE/ffpJ 

where  die  dot  denotes  time  differentiation,  say  i  s  ds/dt.  The  dynamical  metric  has  die  matrix 
dements 

d^lnS/dz*dz^  =  (IL7) 

and  die  total  energy  E  is  a  sum  of  die  electronic  energy  (ijEx)/(s|s)  {fi  being  the  dectrmnc 
Hamiltonian),  the  nuclear  kinetic  energy,  and  the  nuclear-nuclear  rqn^on  energy. 

The  intention  this  very  large  set  of  coupled  differential  equatirms  in  time  is  a  technical 
ptddem  diat  has  been  sdved  in  practice  for  finite  molecular  systems.  The  actual  dioke  of  baris 
functions  and  electronic  preameiers  a  are  discussed  in  sections  V  and  VIL 
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HI.  Sptclai  eonsicterationt  for  extendod  systems 


Long  end  thort  range  potentials 

bi  Older  to  treat  an  extended  ayatem  in  a  nreaningful  way  one  cannot  aimply  let  the  number 
ci  nuclei  and  electrona  grow  indefinitely.  Great  care  ia  required  in  order  to  insure  that  the 
expeeaaionaacaleconectly  with  aize  and  diat  electric  neutrality  is  maintained.  This  is  particularly 
important  anoe  die  Coultanb  interacticm  ia  l<mg  range.  The  6aal  quantities  are  always  mq;nessed 
per  unit  vtdume,  per  unit  cell,  per  atom,  or  per  electron  and  thus  pn^ter  scaling  with  size  is 
crudaL  This  most  be  kept  in  miud  udien  partitioning  die  energy  or  other  quantities  into  partial 
quantities.  Eadi  part  must  be  proportional  to  M,  die  number  of  unit  cells. 

It  is  well-known  (aee  e.g.  Fetter  and  Wakscka  [Ifi])  dial  transiticm  firom  a  finite  to  an 
extended  system  with  Coulmub  interactums  ia  extremely  sensitive.  It  may  be  handled  by  means 
oi  two  coiqiled  limiting  procedures.  One  rqilaoes  the  true  Coulmnb  intoaction  by  a  scremied 
inteiactioo 

(HLl) 

When  the  derired  quantities  have  been  calculated  (»e  lets  the  screening  constant  in  the  cmie- 
qKmding  eiqweasions  tend  to  zero  at  the  same  time  as  the  aze  (the  number  of  unit  cells)  of 
the  system  tends  to  infinity  in  such  a  way  that  the  interaction  remains  operative  throughout  the 
whole  system. 

Requirements  of  diis  type  imposed  tty  die  {diysics  of  the  {woblem  also  show  up  in  the  form 
ci  die  need  for  lattice  sums  to  be  convergent  partititming  of  the  total  enmgy  into  three  parts 
(qiatt  fitom  die  nuclear  kinetic  energy,  ^primdi  by  itself  is  proportional  to  the  total  number  of 
nuclei)  is  an  important  example;  each  part  chosen  proportional  to  the  total  number  of  electrons 
in  the  system  [2S].  Thus,  E  ^  Ei E2  •¥  Ei,  uHtiere 

'  -5/ (ni-2) 

is  the  electiooic  kinetic  energy. 


7(616) 


2  J  ri2 


(in.3) 


is  die  ram  tti’  die  nuclear-nuclear  rqnilsion  energy,  die  etoctrcm-nuclear  attracticHi  enmgy,  and 
the  dectron-electrai  Coolmnb  r^ulsion  energy,  Le.  the  total  electrostatic  energy^  and 


d^id^s 


(IIL4) 


is  the  exchange-eonelation  energy.  In  these  expressimis  7  is  the  first  order  and  T  the  second 
order  reduced  density  matrix  d  the  electronic  system  in  the  state  under  consideratitm,  ^  is  the 
duttge  ci  nucleus  (core)  g,  and  rig,  and  are  die  iqiprofniate  interparticle  distances. 
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Ite  dectricd^  neutnl  system  has  M  nuclei  and  N  elections  so  that 

The  one  matrix  is  pariitiocied  such  that 

-rfeKl)  -  ETtCeiKl)  a>i-« 

i-1 

Le.  each  term  is  associated  with  one  atom  (one  unit  cell).  The  chaise  associated  with  each 
term  is 

imi^ing  diat 

M 

N  *  ^  ng.  (IIL8) 

#■1 


Eadi  one  of  the  three  teems  in  die  dectiostatic  eneisy  (in.3)  contains  a  Imig  lange  component, 

which  is  such  that  the  three  contributions  canoeL  The  fiist  term,  the  nuclear  repulsion  energy, 
is  purely  long  range  and  will  be  left  as  it  stands.  In  order  to  separate  the  long  and  short  range 
contributions  to  the  odier  two  die  potential  doe  to  die  diaige  distribution  jgi 

v,(n,) = / 

is  analyzed.  This  is  a  funedon  of  ri.  Eiqiaiiding  l/r^  in  spherical  harmonics  and  reairanging 
one  obtains 

(n#)  ®  ^  “  w(nf )»  (OLIO) 


triiere 


(n#)  =  y»’2fd»-3#  J  dOa,  J  dCa7#U2l6)l;^  - —1 


dlLll) 


Here  r<  is  die  smaller  and  r>  the  greater  of  rif  and  ra,,  reqiectivdy,  dia  the  angle  between 
n,  rad  t^.  and  dfla^  die  volume  element  relating  to  die  polar  ra^  of  r%. 
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Ute  of  00.10)  leads  to 


fyiS  t 


Sinular  manipoladons  of  die  electron-niickar  attncdcm  term  ^ves 


7>(6I6) 

>7f 


=  -  E  -  E  ^»  /  <«• 


(10.12) 


00.13) 


Therefore  die  long  range  ocnnponent  of  the  electrostatic  part  £lz  of  die  total  energy  can  be 
wiitten  as 


00.14) 


which  obviously  vanishes  when  ^  «  iig  for  all  g. 

The  short  range  part  of  the  electrostatic  intetacdmi  can  then  be  eiqitessed  as 

+  5EK*^»-"*M4,»)-  f  <«i7,«iI6V«»)1- 

nt* 


00.15) 


Hartree-Fock  ▼emu  density  ftmctional  theory 

Whenever  extended  systems  are  dtscossed  it  is  unavoidable  to  touch  the  question  of  die 
bask  ooe-electton  ^iproximadmi.  It  becanre  dear  very  early  that  the  Hartree-Fock  method  is 
unsuitable  for  large  systems,  particularly  for  nmtals  whoe  it  leads  to  a  vanishing  dmisity  of  states 
at  die  Fermi  energy  [26].  This  is  easily  semi  for  the  electron  gas  and  it  has  been  shown  much 
later  bow  this  pedidogy  is  the  result  die  ctmjuncticm  of  three  factors  (i)  Coulombic  forces,  (ii) 
extended  systems,  and  Oil)  die  restricted  Hart^Fodc  method  [27,  28].  For  nearly  thirty  years 
dendty  fhocdonal  mediods  [29],  baaed  on  the  Hdienberg-Kdm  thecnem  [30]  have  provided 
a  beam  zerodi  order  description  of  extended  systems,  whidi  has  led  to  an  unprecedmited  and 
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extremdy  useful  interaction  between  theory  and  experiment  In  the  practical  implementations 
of  densiQr  functional  methods  the  so  call^  ‘local  dmisity  ^>proximation"  of  ^  exchange- 
cocielatirm  prMential  still  dominates.  Hie  weakness  of  that  approximation  together  with  the 
difficolQr  of  making  systematic  improvmnents  have  during  the  last  few  years  lead  to  somewhat 
of  a  renaissance  of  the  Hartree-Fock  method  for  extended  systems  (see  for  instance  volume  42, 
issue  1  (1S)92)  of  the  International  Journal  of  Quantum  Chemistry).  This  renewal  is  coupled  with 
the  awareness  of  the  necessity  to  combine  Hartiee-Fock  with  corrections  to  remedy  pathologies 
as  the  one  alluded  to  here.  In  the  electron  gas  case  it  has  been  shown  explicitly  how  the  Random 
Phase  Apinoximation  (RPA)  can  provide  such  corrections  [31].  This  is  sometimes  summed  up 
by  saying  that  *llarttee  is  a  better  starting  point  than  Hartree-Fod^*  for  extended  systems  [32]. 

It  should  be  noted  that  although  the  END  approach  employed  here  uses  a  single  determinantal 
wavefunction  for  the  electrons,  it  is  not  synonymous  with  the  Hartiee-Fock  method.  The  time- 
dependent  END  theory  at  this  level  of  sfiproximation  can  be  labeled  as  being  fully  nonlinear 
Hnre-Depeiident  Hartree-Fock  (TDHF).  The  RPA  is  equivalent  to  linearized  TDHF  and  thus 
the  present  description  contains  the  necessary  corrections  to  the  pathologies  observed  for  the 
Hamee-Fock  method. 

IV.  The  model 

The  purpose  of  the  present  paper  is  to  apply  the  END  thetny  to  an  extended  system.  This 
requires  a  model  that  is  sufficiently  realistic  yet  not  too  complicated.  A  Hamiltonian  inspired 
by  the  PPP  model  seems  to  meet  those  requirements.  F(v  ^plications  to  extended  systems  an 
important  advantage  of  the  PPP  modd  is  that  it  correctly  describes  the  long  range  aspects  of 
the  (Coulomb  interaction  [33]. 

Originally  the  PPP  method  was  devised  to  describe  the  (optical)  r  electrons  in  planar  organic 
molecuks.  The  model  employs  a  basis  of  atomic  pz  orbit^,  one  on  each  carbon  atom.  In 
order  to  simplify  the  calculations  one  invokes  the  so  called  differential  overlap”  (2IX)) 
ap[m2ximation,  meaning  that  in  the  many-center  integrals  charge  densities  associated  with  atomic 
orbitals  on  two  different  atoms  are  neglected.  This  can  be  considered  a  reasonable  iqiproximation 
if  the  atomic  mbitals  (AO)  ate  interpreted  as  “orthogonalized  atomic  orbitals”  (OAO)  [25,  34]. 
On  the  other  hand  OAO’s  can  be  regarded  as  Wannier  functirms  [25]. 

As  in  the  PPP  model  the  Hamiltonian  is  expressed  in  second  quantization  and  associated  urith 
a  certain  set  of  basis  spin  orbitals.  For  present  purposes  this  basis  can  be  the  (exact)  solutions 
to  the  RHF  problem  at  the  equilibrium  positions  of  the  nucleL  The  canonical  solutions  ate  then 
Bloch  functions,  whereas  the  cmiesponding  Wannier  functions  are  certain  unitary  transformations 
of  these  Bloch  functions.  The  corresponding  LCAO  (linear  combination  of  atomic  orbitals) 
approximation  would  employ  linear  combinations  of  OAO's  for  the  Wannier  functions  and  linear 
combinadmis  of  Bloch  sums  of  the  basic  AO’s  for  the  canonical  solutions.  In  eitho:  case  one 
teficts  to  the  exact  or  approximate  solutions  of  the  RHF  problem  at  the  equilibrium  nuclear 
positicHis  as  SCF  Bloch  fimctimis  and  SCT'  Wannier  functions.  Tbe  reasons  for  this  choice  of 
basis  are  discussed  in  more  detail  in  section  V. 

The  SCF  Wormier  functions  (WF)  ate  expressed  as 

Mm4r(0  =  u(F-  »5»MC)  =  «m(»0<^(C)  (IV.l) 
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(IV.2) 


inchiding  a  tfia  factor,  and  die  SCF  Bloch  fimctions  (BF)  bectmie 

M/2-l 

».(*,{) =v(t,o»«)  E 

mmii/2 

The  diverse  tiansfcMnnaticm  (i^  *  linc/Ma)  can  be  expressed  as 

M/2-l 

=  E  av.3) 

KmJi/2 

Ibis  means  diat  die  wave  number  Jk  lies  in  die  interval  — x/a  <  k  <  x/a,  which  is  the  first 
BiiUouin  acme  (BZ).  The  electron  fidd  operatm  {a,(m),at(m)}  ctmespond  to  the  WF  and 
{kvik)^  ^(1;)}  to  the  BF,  where  die  subset^  r^srs  to  ibt  Qun.  The  creation  operators  transform 
as  die  ^fin  odntals 


3K 

»*(»)  =  W-'Z* 

m 

BZ 

«l(m)  =  Y,  klik)e-'*^*”^ 

k 

and  the  adjennt  relatkms  hold  for  the  annihilatitm  operators. 

The  Hamiltonian  in  die  WF  rqiresentation  can  be  expressed  as 

BK  BK 

^  *  -Ho  +  53  53  ”)4(”»)o«T(n) 

j  BK 

2  ^  S  »»»(m)n^(n) 

*  m,n 


(IV.4) 


(IV.5) 


where  the  occupation  number  operator  n,(m)  s  at(m)a«r(m)  and  where  in  the  last  sum  the 
term  widi  m  «  n  and  a  ^  vk  missing. 

The  diagmial  matrix  dement  a(m)  of  the  one  Hamilttmian  with  respect  to  die  WanniOT 
function  centered  at  m  is  (Q  is  the  atomic  cme  charge) 


BK 

o(m)  s  a„  - 

mikm 


Q 

|^(m)  -  nl 


(IV.6) 


The  first  term  dqiends  on  the  cme  displacements  p{m)  (t^iich  are  not  restricted  to  any  one 
directimi)  via  a  shent  range  potential  The  sum  can  be  separated  into  two  terms.  A  Coulomtnc 
one  ind^endent  of  die  can  displacements  dud  cancels  amiiar  terms  in  the  dedronic  and  the 
nuclear  rqnilsiaa  energies  and  one  that  depends  on  the  core  displacmnmits.  The  off-diagonal 
dementi  ^(m,n)  of  die  core  Hamiltonian  also  depmid  cm  the  displacements  via  short  range 
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pomtials.  Hie  ptnunelen  j(n  ->  m)  lepiesent  ^  interaction  between  electron  density  on  atom 
m  (or  density  in  cell  m)  tod  t^  on  atom  II.  Due  to  the  choioe  of  basis  functions  these  parameters 
d^end  on  die  distance  |f?i  -  =  A  but  not  cm  the  displacements.  Following  Stolarczyk  et 

aL  [35]  thefonn 

is  diosen  to  obtain 

X-\Jl)  *  -i2(l  + 

which  shows  dial  the  term  A(il)  is  free  of  spurious  Imig  range  effects. 

The  term  J7o  in  the  Hamiltonian  rqiresents  the  rqiulsion  between  cores 

If. -Ly'  0* 

*  1^™)  -  ^(“)l 

As  discussed  in  secdtHi  n  the  total  wavefiinction  is  assumed  to  be  a  single  Slater  detmninant 
The  cmesponding  one  matrix  or  the  Fock-Dirac  density  matrix  in  the  Wannier  lepiesmitatirni 
can  be  expressed  as 

BK 

m,ii  9,0* 

The  matrix  elements  n)  clearly  depend  on  the  choice  of  basis  and  constitute  an  avmage 

over  the  state  under  ctmsideration 

D»»'(m,n)  =  (aj,(n)a,(m)).  (IV.ll) 


(IV.7) 

(IV.8) 

(IV.9) 


The  notaticm 

D*®(m,  n)  +  D^^(m,  n)  =  D(m,  n)  (IV.12) 

is  introduced  demanding  rhat  m)  s=  m)  =  QI2  in  the  refmence  state,  whme  Q 

is  the  core  charge  of  the  only  type  of  atom  presmit  in  the  chain. 

The  average  value  of  the  various  terms  in  the  Hamiltonian  requires  the  evaluation  of  the 
expectation  values  (n,(m)n^(n)},  which  for  a  single  detmninant  becomes  [20] 

(n,(m)n^(T»))  =  (n,(m))(n<r'(n))  =  m)D®'‘^'(n,  n).  (IV.13) 


The  eaqiectatitm  value  of  the  Hamilttmian  (IV.S)  with  respect  to  a  single  detenninantal  state  can 
then  be  eiqnessed  as 


a  bk 


1 


BK 

+  ^  a{m)D{m,  m) 


m 


2 

BK  BK 

+  7(0)  D®®(m,m)D^^(m,m)  +  53  n) 


av.14) 


m 


-  BK  fin 

+ 5  E  T(»  -  "■)  E  >»)^'''(’».  >•)• 
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(IV.15) 


Uaiof  die  above  definitiona  inclodim  (IV.7)  and  die  eaqianaon 


|.5  —  R 

for  r  «  die  numopole  terms  are 


1  1  a? 


2R^ 


.  BK 


[g  -  D{m,  m)l[Q  -  Djn,  n)] 
jm-nj 


(IV.  16) 


This  ejqxesskm  corresponds  to  (lELU)  in  the  general  treatment  ¥ot  ionic  systems  this  is  the 
Madelung  energy.  It  vanishes  when  Q  s  D(mym)  for  all  m,  a  condition  that  means  that  the 
number  of  dectrcms  associated  with  atcm  m  is  die  same  as  the  con  charge.  This  is  the  case 
for  die  reference  state.  When  the  nuclei  vibrme  and  the  electrons  react  to  this  motion  the 
atomic  charges  will  fluctuate.  The  fluctuations  around  the  equilibrium  charge  add  up  to  zno, 
while  oquession  (IV.16)  rqnesents  the  Madelung  energy  of  these  fluctuating  atomic  charges. 
It  fdlows  duU  the  requirement 

BK 

5^  D(m,  m)  *  MQ  (IV.17) 


is  a  reasrmable  one. 

In  sectitms  Vn  and  Vm  a  particular  kind  of  fluctuations  ate  introduced  such  that 

/)(m,m)  =  g  +  (-)”*^, 

which  obviously  satisfies  condition  (IV.17),  and  gives  the  Madelung  energy 

2a  "  |m  — n|  o’ 


rnytn 


(IV.18) 


(IV.19) 


This  ccmdition  is  assumed  to  hold  in  the  following  analysis. 

The  next  set  of  terms  in  the  multipole  expansion  (IV.IS)  are 

BK 

Q  n)p.(m),^  _ ,, 

^  4  ■  W  - 

an  expressitm  dut  also  vanishes  under  the  condition  Q  =  I7(m,m)  for  all  m,  and  also  with  the 
qiecial  dhoice  (IV.18)  this  sum  is  zero,  since 


(IV.20) 


|m  —  n|* 

^  ««  BK 


mfiii 

BK 


(IV.21) 
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However,  a  state  witii  general  fluctuating  charges  will  contribute  to  the  Madelung  energy. 
Sirnilaily.  some  of  the  terms  piopoctkmal  to  vanish,  so  that  for  small  vilMations 


fff\  «  V  PwMpM  -  2p»(m)p,(n) 


(IV.22) 

+  7(0)  53  m)D^^(m,  m)  +  V  n) 

m  nytflt 

j  BK 

+  2  £  “  f>^Wrn,m)D{n,n). 

mftii 

In  smneu^iat  mote  detail  than  expressed  by  (IV.18)  one  finds 

X)±(m,m)  =  2  +  (_i)-4t  (IV.23) 

vdiete  the  notations  D""  s  and  s  D-  have  been  introduced.  One  may  then  also  write 


BK 


BK 


(IV.24) 


The  parameters  and  /9(m,n)  are  doe  to  the  short  range  forces  within  and  between  the  cores 
and  ate  expanded  in  terms  of  the  displacements  as 


Qohn  =  Qa®  +  53  +  2  S 

i  iJ' 


(IV.25) 


and 


^(n,m)  =  )3®(n,m)  +  53 

i 

2  *.f.» 

+2  S  m). 

iJ' 


(IV.26) 


This  result  has  been  obtained  using  the  tetatum  /3^(n,m)  =  n)  for  the  first  derivatives. 

The  total  energy  separates  into  two  parts,  one,  JE^  consisting  of  those  terms  that  dq)aKl  on  the 
nuclear  di^Ucements  and  another,  bdng  die  energy  of  the  system  with  the  static  lattice,  le. 


(IV.27) 
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(IV.28) 


2  M 


«  M{Qa^  +  Y,  »)  +  -rWIx  +  2  S 


M^O 


■®'‘  =  ^  E  ji~J^I»  !<’»("■) + '^(’”) + '^(”*)1 


OmC 

.^’Sr  />«(”»)/>«(”)  +  />#(”»)/>»(») - ^pM{m)p,{n) 

■^2a*^  Im-n|a 

•,»,*  ,  BJf  *^,* 

+  52  ]C  +  2  5^  ^  ^ii*/»iMPi'(»»») 

*•  i  *»  i.i' 

BJT  «if.x  .  c,f,x 

"H**  i  iii' 

(IV.29) 

The  fifst  and  die  last  term  in  die  eaqwcssicm  for  JEtf  are  really  of  higher  Older  and  are  not  part  of  die 
harmonic  aptnofrimatkm.lwit  kept  here  to  riiow  how  such  terms  enter.  The  nockar  di^lacemeots 
and  die  d]fnamical  vaiiablea  z  can  be  chosen  such  diat  da^  are  aero  at  a  stadtmary  point,  vMch 
means  diat 


djEft  -i-  Bj)  j  __  d{Bgt  -f  Ej)  I 

dz  |**0~  Qp  |z  =  0 

ps=0  pssO 


(IV.30) 


are  the  equations  defining  die  stadonary  point  The  symmetrized  factor  [D(m,n)  +  D{n,m)] 
gaaranlBes  a  real  result  for  die  last  term. 


V.  Elactronic  basis  functions 

The  dioioe  of  basis  functions  for  the  description  of  electrons  in  a  systmn  with  moving  nucld 
requires  careful  analysis.  Only  the  situatitm  i^iete  die  nuclm  undogo  small  vibrations  always 
keqiing  diem  dose  to  dieir  equilibrium  positions  is  ctmsidaed.  In  piiiKtiple,  a  complete  set 
of  electnmic  basis  functions  is  needed.  Formally  such  a  set  is  available  in  the  form  of  all  the 
sohitmits  oi  die  band  themy  problem  for  the  equilibriom  positions  of  the  nucleL  Such  a  basis  is 
obtained  as  the  solutions  to  an  effective  one-electron  equatitm  for  a  static,  fiilly  periodic  lattice. 
At  this  stage  cf  die  discussion  there  is  no  need  to  specify  that  equation  in  more  detail.  It  could 
be,  say,  of  Hartree-’Fodt  type,  of  density  functional  type,  or  somedimg  elre.  What  is  needed  is 
die  property  common  to  all  such  equatitms,  namely  ^  translational  symmetry,  which  implies 
that  the  aolutions  are  Blodi  functitms,  Le.  eigenfii^tms  of  the  translation  qierator. 
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In  actual  CialcalatioBi  a  ^peci6c  choice  of  such  Blodi  functuMU  obviously  has  to  be  made. 
The  calculatioa  d  die  enetgy  bands  and  the  cwiespondiiig  wavefiinctimis  is  a  wdl  developed 
proceduie  [36, 37]  and  diefe  are  a  number  of  possible  choioes  for  rqnesenling  the  solutions  to  the 
stttic  lattice  problem  in  terms  of  i^iecific  ba^  fiinctitms.  This  pieseatatum  limits  the  discussion 
to  tte  case  when  die  original  basis  fhncdons  are  atomic  oibitals  of  some  kind  centmed  at  the 
eipiilibriumnnclearposidoiis.  Rom  each  type  of  atomic  cubital  in  a  unit  odl  a  Bloch  sum  can  be 
formed  adapted  to  tte  translational  symmetry.  The  final  solntions  of  die  effective  one-electron 
proUem  are  dieii  linear  ccunbinations  of  such  Bloch  sums.  The  term  *SCF  Bloch  functicms'  are 
used  for  sudi  solutions. 

Given  any  set  of  Bloch  fimcticMis  in  a  band  the  cmresponding  set  of  Wannier  functicms  can 
always  be  constructed  by  means  of  a  unitary  transfcumation  [38].  If  that  is  done  with  a  set  of 
Blodi  functions  eiqiressed  as  linear  combinations  of  Bloch  sums  of  atomic  orbitals,  dm  resulting 
Wannier  functions  can  be  closely  approodmated  by  linear  ccunbrnaticms  of  ortitogcmalized  atcunic 
orbitals.  But  die  unitary  transformation  can  also  be  carried  out  on  the  exact  solutions  to  the 
static  lattice  problem,  ti^  may  be  valuable  from  a  formal  point  of  view.  The  formally 
exact  Wannier  functions  are  then  obtained,  uhirii  ccmstitute  a  set  of  non-canonical  mbitals  for 
die  single  d^ermiiumtal  wavefiinction  to  rquesent  the  solution  of  die  static  lattice  problem. 

The  functicms  (IV.I)  and  (IV.2)  {uimatily  rqnesent  sudi  exact  canonical  and  ncm-cancmical 
sdlnticms  of  the  effective  one-electrcm  static  lattice  {Moblem.  They  do  not  ccmstitute  a  cmnplete 
set  of  functions  sinoe  they  correspond  to  <mly  (me  band.  But,  in  principle,  it  is  possibte  to  include 
more  dian  one  band.  These  functions  provide  a  good  descr^on  of  the  syston  at  die  equilibrium 
nudear  positions.  For  smsU  departures  fitmn  die  equilibrinm  positions  it  would  therefore  semn 
natural  to  emid(ty  dieae  fimctim  as  a  basis  for  die  END  treatment  Since  die  term  *hasis 
function"  is  used  in  several  different  c(mtexts  one  should  clarify  the  following  points.  The 
dectroinc  basis  fimctions  for  die  END  treatment  are  the  exact  m  i^iproximate  soluticms  of  the 
effective  one-dectrcm  equaticm  for  the  static  lattice  problem,  Le.  functions  (IV.I)  and  (IV.2).  The 
parametma  diaracteiizing  die  FFP  Hamiltonian  are  thus  expressed  in  terms  of  these  functicms. 

hi  ab  initio  cakulaticms  obviously  mote  specific  dioices  have  to  be  made.  One  can  discon 
two  dittinct  ahematives.  One  would  be  to  use  die  primitive  basis  of  Gaussians  centered  at  die 
nuclei  These  basis  functions  ate  allowed  to  fcdlow  die  nuclear  motion.  At  any  given  moment 
die  nuclear  firamewoik  shows  no  particulv  symmetry.  This  apiuoach  has  led  to  succesriiil 
descrqitioos  of  small  molecular  systems.  However,  for  an  extended  system  die  lack  of  periodkity 
means  diat  one  cannot  speak  about  Blodi  or  Wannier  functicms,  leading  to  great  complicaticms. 
The  (Mher  possibility,  uhich  has  been  chosen  here,  is  based  cm  the  notion  diat  the  nuctei  never 
dqiart  very  far  firom  die  regular  equilibrium  positions  and  the  SCF  Bloch  or  Wannier  functicms 
will  be  adequate  dioices. 

Using  the  Bloch  or  Wannier  functicms  in  one  band  means  that  the  possibility  for  die  dectrcms 
to  adjust  to  die  changing  nuclear  envinmment  resides  in  mixing  the  qnn  orbitals  used  in  the 
reference  determiiuttit  with,  in  principle,  all  virtual  spin  orbitals.  In  ordn  to  eiqilore  the  effects 
d  sudi  mixing  in  a  scmiewfaat  systematic  fashion  die  pairing  of  eadi  reference  qnn  orbital  with 
one  virtual  qrin  oritital  is  first  considered.  This  can  be  looked  iqxm  as  a  successive  lowering  of 
die  translational  symmetry.  Since  die  transhukmal  symmetry  can  also  be  used  to  classify  normal 
modes  d  vibraticm,  such  pairings  should  [ffovide  a  number  of  useful  qiecial  cases  of  END. 
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VI.  Lattlct  dynamics 

In  Older  to  connect  die  END  treatment  of  vil»atic»al  and  electronic  structure  with  more 
traditional  procedures  this  section  treats  die  qiecial  case  n^iere  the  total  eaagy  in  the  equations 
of  motion  (ILd)  does  not  dqiend  <m  the  electronic  dynamic  variaUes. 

Widi  oidy  nuclear  ttynamic  variables  the  equadtms  of  modem  reduce  to 

-dP/d»-dE/9R 

dR/(U  =  aE/dP  '  ' 


ie.  die  classical  Hamilton’s  equadons  of  meMion. 

Restricting  die  cemsideradoos  to  longitudinal  vitadions,  Le.  dynamics  for  which  px{m)  = 
Py(m)  s:  0  for  all  m  leads  to  simpler  and  mene  transparent  exixessions  than  those  of  the  full 
treatment  bridle  dearly  displaying  the  capabilities  of  the  model  Adding  the  nuclear  kinetic 
energy  to  (IV^,  and  using  (IV.30)  yiel^ 


BK 


E  ss  Egt  +  Ei  +  ^  ^ 


m 


JL 

2itf, 


m^% 


m 


2Mn 


The  core  charge  issettofi’^ltodthe  fdlowing  notation  is  introduced 


(VL2) 


(VL3) 


widi  s  «,y,or  z]  It  should  be  noted  that  die  density  matrix  of  the  reference  state  satisfies  the 
rdatkm  l>o(u,m)  s  I)o(n  —  m,0).  In  equation  (VL3)  and  in  the  expression  fca  Egt  the  drasity 
matrix  Do  of  the  refiBieaoe  state  ^qiears.  The  energy  derivatives  with  respect  to  displacmnents 
dmi  are 


dE 

dp,{m) 


(VL4) 


Thus,  die  coiQiliitg  between  different  nuclear  displacements  is  due  to  the  ‘Tiopping  parameters" 
jS  in  die  Hamiltoaian  and  also  introduced  via  dm  multipole  expansion. 

The  equations  ci  motion  for  diis  qiecial  case  (longitudinal  nuclear  motion  and  one  land  of 
miclei  widi  mass  M^) 


dPgjm)  BE 

dt  dpgitn) 

dpg(m)  __  dE 
dt  dPtim) 


(VL5) 
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be  oombiiied  into 


(VL6) 


d^ptim)  BE 

"  dti  dp,(my 


Expending  die  di^lacements  in  nocmal  coocdinates  Q(q,t)t  sndi  that 


1 

/»*(*«)  = 


(VL7) 


sabsdtnting  due  into  the  equatirm  of  motkm,  multiplying  by  exp(~t9am),  and  summing  over 
m  yield 


4\  ^  ga>Eii 

-  EB„(n,0)«‘*“W(,,«).  (V18) 


The  Ofuatz 


(VL9) 


used  in  equadcm  (VL8)  yields  the  following  ex{Nression  for  die  Icmgitudinal  frequency 

,  „  BK  BK 

*  '  *  ||a|M 

,  .  jr/a-i  .  .  jr/3-i 

“7r{-4  E  ^=^+>*..+2  E  B..(".o)co.(,.»)). 

"  Mvl  nsl 


(VLIO) 


Li  Older  to  get  the  correct  behavior  for  small  q,  Le.  u{q)  0  die  cemdition 

v/2-1  jr/3-i 

E  ^  +  +  *  E  A.(n,0)  =  0 

lucl  asl 


is  needed. 

The  inverse  transformation 


'■o-m 


(VLll) 


(VL12) 


of  the  diqilaoements  together  with  the  condition  that  pg{m)  are  real,  gives  die  relation  Q*{qyt)  = 
Q(— 9,<).  The  general  mqiression  for  the  displacements  can  then  be  written  as 

P*{rn)  =  +  53(Qi(q,f)coeqom  -  Q2(9,t)  sin  qam]},  (VL13) 
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I 


with  the  real  and  imaginary  parts  of  the  normal  coordinates  given  by 

^i(9i  <)  =  Qoi(9)  +  ^02(9)  nnun{q)t 

Q2{q,i)  *  -Qoi(9)sinwi(g)t  +  Qn{q)coBun(q)it. 

Should  transversal  vitodioos  also  be  included  in  (VL2)  die  general  oq^ession  becomes 

BK 

E^E  +  ^  ~ 2p,(m)/>,(n)  ^ 

**  2e?  "  Im  ~  nl* 

Mftm  ' 

+0  ]C  S 


(VL14) 


(VL15) 


•"  iJ' 


BK  «,f,* 


«  ••n  '»»»«  BK  02 

+2  S  ^  + 23  2]ir' 

Mflm  iJ'  m  ■ 

Urns,  die  one-oenter  parameters  induce  couplings  between  different  nuclear  displacemmit  com¬ 
ponents  on  the  same  atom,  u^iile  the  'liopj^  parameters**  0  produce  the  couplings  between 
all  compcments  on  different  atoms. 

VII.  TIme-dapendent  band  thaory 

As  anodier  ^wcial  case  die  END  equations  for  fixed  nuclm  ate  studied.  This  is  the  Bom- 
Oppenheimer  starting  point  for  most  problems  of  dectronic  structure.  The  equations  of  motion 
(IL6)  dud  determine  tte  time  evolution  of  die  electronic  state  now  reduce  to 


fiC  0  \f  dx/d<  \  _  fdE/dz*  \ 
\  0  -iC* ) \iz*/dt dE/dz  )' 


(mi) 


The  electronic  basis  crmsists  of  all  the  SCF  Bloch  (IV.2)  or  SCF  Wannier  (IV.l)  functions  in  a 
particular  band.  The  most  general  combinatirm  of  qnn  cubitals  is  of  the  form 

RZ 

•>.(*.£) +  (ttoLZ).  (vn.2) 

Here  the  abtueviation  LZ  s  Little  Zone  is  used  for  die  interval  -ir/2a  <  k  <  ir/2a  and  RZ 
s  Residoal  Zone  for  the  remainder  of  the  first  Brillouin  Zmie  (BZ),  which  consists  of  the  two 
intervals  —v/a  <  k  <  — ir/2a  and  irf2a  <  k  <  ir/a.  The  r^»mice  determinant  is  doubly 
filled  widi  the  orbitals  v{ktf)  in  the  LZ.  Ihe  virtual  ortntals  v(Jk',f)  are  charactnized  by  the 
wave  numbers  if  m  RZ. 

The  possilHlity  of  using  fully  general  spin  cnbitals  like  (V1I.2)  should  certainly  be  considered. 
Here  die  case  wheae  each  spin  orbital  has  dther  am  0  spin  is  first  eiqilrned.  The  corresponding 
orbitals  are  then  in  gaiersi  difident  Qim  set  of  orbitals  w^.{k^i^  is  combined  with  a  and 
another  set  ti;_(ib,f)  with  0  spin, 

RZ 

jiz 

U7_.(*,f)*»(*,f)  + (*  in  LZ). 

V 


16 


Febroary  23. 1994 


Eadi  tet  liieae  Amctioot  can  widiout  loss  of  generality  be  cboaen  otthogwal,  even  thoagb 
Aby  are  not  nocmaliwxi  to  unity,  le. 


J  w^{k,  r^vt±{U  rf)d?  *  d±*^« 

4±*  «  1  +  *±,k>k*±;k»k  =  1  + 


(Vn.4) 


Becauae  of  die  pairing  dieofem  [^]  die  z  parameters  can  be  assumed  to  satisfy 


J  w%{k,r)u}:f{l,?)df^X±kSu 


(VIL5) 


Thus,  both  die  two  (Iif2)x(i02)  matrices  and  and  die  product  matrices  z^z_  and  zLz+ 
are  diagooaL 

Ibe  coneqionding  Fock>Dirac  density  matrix 


k 


(vn.6) 


can  dien  be  expressed  in  terms  of  die  z*s  (see  e.g.  [Q).  The  basic  Bloch  functitnis  are  here 
assumed  to  be  aotutkxis  d  a  restricted  Hartree*Fodc  problem  and  time  reversal  symmetry  would 
imply  diat 

^-k,f)  -  v‘{k,F) 

There  is  no  reason  to  impose  such  relaticms  on  the  functicms  (VIL3). 

In  order  to  iUostrate  how  the  END  method  works  die  case  widi  only  two  z-parameters  per 
wave  number  in  LZ  is  studied,  so  diat 

w+(fc,  f)  =  v{k,  f)  +  z+.jk.fcw(fc',  f) 
w_(fc,rO  =  v(fc,f)  +  z^.^kv(k\r). 

This  is  an  example  d  pairing  whidi  can  be  made  in  several  different  ways,  both  with  respect  to 
die  qiatial  orlntals  and  widi  reqiect  to  qiin  [40].  A  determinant  with  half  of  the  electrons  filling 
w+  orititals  widi  qnn  a.  and  the  <xber  half  filling  v>-  orbitals  with  si»n  forms  die  basis  for  the 
Alternant  Molecular  Orbital  (AMO)  method  (see  [41]  and  refismces  therein).  In  applications  of 
the  AMO  mediod  to  the  description  of  solids  [42]  the  pairing  is  chosen  sudi  diat 


(Vn.8) 


y^k  +  rla;  ib  in  LZ;  in  RZ. 


(Vn.9) 


One  could,  of  course,  leave  V  unqiecified  within  RZ  in  developing  the  formal  theory,  and  a 
vati^  of  dioioes  could  be  eaqdored  in  actual  applications.  In  order  to  draumstrale  die  principles 


17 


FBtaiUfy23. 1994 


of  die  END  dieofy  it  sufBoes  to  limit  the  pieaem  development  to  the  choke  in  equati<Hi  (V1I.9}, 
«Bd  write 


*  *±k- 


(mio) 


This  meutt  that  there  are  W1  dynamical  variables  a4k(t)  for  ekctrons  with  op-spin  and  an 
equal  number  of  variables  «-s(t)  for  down-si^  electrons.  The  density  matrix  coefficients  (le. 
dements  of  the  so  called  charge  and  bond  order  matrix)  are  dien 


sudi  that 


k 

+ (-!)"»*»  +  K-i)”-  -  11I»±*P  +  •••}. 


D(m,n)  ss!  D^(m,n)  +  P_(m,n)  = 

n)  -I-  n)  -h  I>2(m,  n) 


(mil) 


(mi2) 


in  terms  of  the  zeroth,  first  and  second  mder  terms  (terminating  after  secrmd  order)  in  the 
dectronk  ^namkal  variaUes.  This  gives 


BK 

5^  ^(m,n)D(m,n) 
»  »K 


_  .  «  ^  ■ 


(mi3) 


A 

E 


-l)/2 


-/S*’(m,0). 


m»odd>0 


For  the  diagonal  elemoits  in  (VILll) 


iJ±(m,m)  =  i  +  ^{4,  +  »«}  -  i  +  (-!)• 


(mi4) 


(mi5) 


Adding  die  contributions  from  the  two  qrins  yields 


D(m,m)  =  D^{mym)  -I-  D^im^m)  =  1  +  (—1)"*? 


(mi6) 


widi  +  aid  with  (-l)”'9(t)  rqnesenting  die  fluctuating  charge  at  site  m,  such  that 


BK 

£(-i)-,W  =  o 


(mi7) 
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at  ail  timea.  An  eqnalkm  for  the  total  atatic  energy  form  (IV.28)  can  then  be  obtained  as 

InO  1 

V  +  5  £  (-» +  “• 


BK 

+  B  +  7(0)It  +  ?+«-)  +  5  ^(”)) 

BK  IZ 

-2  £  + 
maodd  k 


DiffiBeentiation  with  respect  to  the  z  vaiiables  gives 


^  =  2A)«  +  7(0)4-  -  2B(k)z%i 
^  -  !M.4  +  7(0)4+  -  iB{k)zU 


and  dieir  conqdex  conjagates,  s^iere. 


(VIL18) 


(VIL19) 


(V1I.20) 


B{k)^  ^  i9®(m.0)c“*"». 


(VE.21) 


Fbr  kinLZii  holds  diat  B{k)  <  0.  It  follows  fitom  (V1L19)  that  the  derivatives  vanish  when 
the  t  variables  do,  Le.  the  electronic  dynamical  variaUes  are  defined  relative  to  a  stationary 
point  of  Egi, 

The  dynamical  metric  is  obtained  fitom 


giving 


s  -  {.w  -  n «» + 

k€LZ 


(VIL22) 


^  S^hiS  c  u.  .2^2  (vn.23) 

C+fc,_|  =  C'_*,+|  as  0. 

For  the  harmonic  ^>inoximation  the  dynamical  metric  is  needed  only  through  zeroth  order,  Le. 


C+k,-¥t «  C^k,-i  *  Su 


(VIL24) 
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1  I 


yicMfaig  die  equatkns  of  motkm 

=  2M  + 1(0)9-  -  2B(t)»+4 

-  2i4,4  + 1(0)9+  -  2B(k)z.t 

and  thdr  ooopia  conjugates,  bmoducing  die  variaMes  Uk  —  z^k  -*  leads  to 


and  its  complex  conjugate.  In  matrix  torn 


ti  0  ir ui  r 

.0  -ilJlu*J“[ 


Ai-2B  Ai 


Ai  Ai 


ii  lf«l 

-2BJ[uV 


(VIL25) 


(VIL26) 


(VIL27) 


U*  tifc 


a  column  matrix  containing  the  MU  variables  14,  the  (M2)x(M2)  matrix 


Ai  »  - 


7(0)  1  1 


(VII.28) 


(Vn.29) 


and  die  diagmial  matrix 


^(^-jf/4)  0 

0  ^(*-JI//4+l) 


[0  0  ••• 

The  scdudcm  of  the  dynamical  equations  can  be  accomplished  by  the  ansatz 


’ul_rx  life®*  0  Ifvl 

u*J”[l  Xj[  0  e-“**J[v*J’ 


(VIL30) 


(VIL31) 


where  the  time  depmidence  is  confined  to  the  matrix  blocks  which  are  diagonal  with 
elements  exp[±iu;i(ib)t].  Substituting  (Vn.31)  in  (VJL27)  and  writing 


V  =  V*  = 


lead  to  a  set  of  homogeneous  linear  equations 


l: 


(vn.32) 


(VIL33) 
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-•  n  * 


a  =  -xn  -  (Ai  -  2B)X  -  Ai 
b*n-(Ai~2B)-AiX. 

Partiticning  of  (V1L33)  yidds 

[a-ba-‘blV*0 

and  die  fimpieiides  (tnik)  a>e  obtained  finm  the  secular  equation 

det{a  —  ba“*b}  »  0. 


(VIL34) 

(VIL35) 

(VIL36) 


Tbe  dynamical  variaUes  can  then  be  eiquessed  as 

LZ 

(VIL37) 

I 

Equation  (Vn.25)  yields  a  completely  analogcms  treatment  for  gu  —  z+k  +  as  that  just 
cmnpleied  for  ii^.  One  gets 

LZ 

gkit)  =  +  Yi  Yuwie^^*^*.  (VIL38) 

I 

Tbe  original  dynamical  parameters  are  obtained  as 


*-*  =  «*). 


(VIL39) 


Vm.  Combined  electronic-nuclear  dynamics 

Ibis  section  presents  a  simple  case  where  both  dectronic  and  nuclear  dynamical  variables 
are  allowed  to  vary  simultaneously.  Tbe  starting  point  is  die  harmonic  approximation  of  the 
energy  for  longitudinal  vibrations  and  for  paired  ttynamical  qnn  cxbitals  as  discussed  in  die 
inevious  secdtm.  Tbe  core  charge  is  set  to  unity  iQ  s  l),  and  tte  total  energy  of  the  system  is 

E^E^  +  Ei  +  Tn  (VnLl) 


widi 


+  i  ^  (-l)’A(n)l  +  o” 

BK 

+  «+7(o)(j+*+«-1  +  5E'*W) 

iifio 

BK  LZ 

-2  Y,  ^(m,0)£(U+.l»  +  |r-»|»le“*". 

msodd  k 
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»  ^  X  /^(m)  fiMim)pg{n) 


gX  ^  BK 

+  S  +  2  ^ 

in  m 

BK 

+  J3{P("»i»»)  +  i>(n,m)l(/S;(m,nV,(m)  +  ~0i,(m,n)p,(m)pg(n)]}. 


(VnL3) 


^  P*(m) 


2M.  * 


(VIIL4) 


The  deader  matrix  D±{rn,  n)  contains  terms  of  zerodi,  first,  and  second  order  in  the  electronic 
dynamical  variables  [cf  (VILll)], 


Z?±(m,n)  =  Z?o(m,n)  +  Z?±,i(m,n)  +  Z>±,2(m,n). 


(VnL5) 


Since  9  is  of  first  order  in  the  electronic  coordinates  the  first  term  in  (VIIL3)  is  really 

of  third  (xder  overall  and  should  not  enter  in  the  harmcmic  aiyproximation.  This  term  is  neglected 
as  well  as  other  tmns  of  higher  order  than  two  in  the  dynamical  variables.  Differoitiation  with 
re^Mct  to  the  nuclear  displacements  yields  the  ftdlowing  expression 


2  ^  />«(») 


^  =  >45-  +  .4,  +  Aggpgirn) 

»,(m)  a*  |n  -  m|*  ' 


+  po(m,n)  +  f}i(m,n)  +  £)o(n,m)  +  I?i(n,m)]j9®(m,n)  (V1IL6) 

n^m 

BK 

+  E  ^o{m,n)/S®,(m,n)^,{n). 
nfhn 

Assuming,  as  is  dtme  here,  that  the  refeteiKe  state  corresponds  to  a  stationary  point  of  the 
energy  implms  that 

dE 

for  all  m.  The  only  term  in  (VULfi)  which  contains  the  electronic  dynamical  variables  is 


^  n)  +  Di{n,  m)]/8®(m,  n)  = 

i»#m 

+^*W(4*+*i*)}. 


(VIIL8) 
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where 


$(k)  =  A(t)  +  iAW, 
A(t)=  E  ^(».0)co.(*an), 

53  i8j(n,0)  sin  (fcan). 

IMOdd 


(VIIL9) 


Since  the  derivatives  /9|°(n,0)  are  odd  functions  n  ^e(fc)  =  0.  Using  (VIL19)  and  (VIILl)  the 
derivatives  with  xespoct  to  the  z  variables  bectmie 


^  =  2A),  +  7(0)«-  -  2B{t)4t  +  2m  jf  E  (-irz-'W. 

’  #W 

BK 

^  =  2Am  +  7(0), +  -  2B(i)»l4  +  2m^  E 

Wl 


(vm.10) 


and  their  omiplex  conjugates.  These  derivatives  vanish  whoi  all  the  dynamical  variables  are 
zero,  Le.  the  dynamical  variables  are  defined  relative  to  a  staticmary  point  of  the  eattgy.  Rnally, 
die  derivatives  with  respect  to  the  nudear  momenta  are 


dE  P,{m) 
dP,{m)^  Mn  ' 

The  four  coupled  equations  of  motion  are 

^  dE 

*  dt  *  dt  ”  dz±it* 

d/»,(m)  dE  dP,(m)  dE 
df  *"aP,(m)‘’  “  dt  ^  dp,{my 


(VnLll) 


(vnii2) 


Introducing  the  normal  coordinates  (VL7)  the  latter  two  equations  become 

dQ(k,t)  P{k,t) 
dt  Mn 


and 


with 


^  e(k)Qik,t) 


Thus,  the  forces  (VIILIO)  becmne 

^  =  2^,  +  7(0),-  -  2B(lb)<,  +  2^/S.(i)(^(,/o,i)/^/j^^: 

^  =  IMo,  +  7(0),+  -  2B(k)zU  +  2iA(t)(?(T/<i,l)/VSM: 


(VIIL13) 


(VIIL14) 


(VIIL15) 


(VHLlfi) 
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and  didr  complex  conjugates.  The  equadons  of  motion  then  beccnne 

=  2AH  +  7(0)«-  -  SB(i)»+»  -  2iA(t)0(ir/(i,i)/v/^ 

(vin.n) 

=  2A.,  +  7(0)«+  -  -  7iff,(k)Q(rla,t)ly/uK 

and  their  cmnplex  conjugates.  In  this  particular  model  only  the  normal  coordinate  Q{vfa,t) 
couples  to  the  dectronic  motion.  The  other  normal  coordinates  satisfy 

!H.^^^=e(k}(}{k,t),  tjtir/o,  (Vin.18) 

with  the  soluticms  QiKi)  =  with  die  firequmicies  givmi  by 

u\k)  =  -e{k)IMn.  (VIIL19) 


niis  is  the  same  ejqnression  as  for  die  longitudinal  frequencies  in  (VLIO).  The  longitudiiud 
frequency  in  (VIEL19)  can  be  expressed  as 

J//2— 1  M/2—1 

£  ‘^1**’"*+^“  +  ^  E  ^o(0.n)^;.(0,n)c«(t«n)].  (VIII.20) 

•  IWBl  dsl 

Correct  bdiavim-  at  small  Le.  for  u>{k)  0  requires  that 


(Vra.21) 


Comparison  of  (VIL25)  with  (VIIL17)  shows  that  the  difference  variable  —  z^k  —  ^-k 
sdll  satisfies  (V1L26)  as  in  the  case  with  frozen  nuclei  Thus,  obviously 


LZ 

z+jk  -  z_jb  =  +  ^Xkivie^^^*, 

I 

as  in  (VIL37).  Introducing  the  abbreviated  notations 

Q  =  Q(t)^Q(TMIy/UM: 

P  =  P(t)  =  F(,/c,<)/VAf^ 
Straightforward  algebra  using  (VIIL17)  gives 

=  (M,  +  7(0))?  -  2B(i)»  -  «A{ife)(3 


(vra.22) 


(VIIL23) 


(Vin.24) 


for  die  sum  variable  gu  —  +  z-a.  Tlie  definition  of  q  in  equations  (VILIS)  and  (Vn.16) 

leads  to  die  eqxessitm 
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with  A  K  (8i4o+27(0))/Af .  This  means  diat  the  leal,  grk,  uid  imaginary,  gut,  parts  of  gt  satisfy 


grk^-2B(k)gn-Afi,(k)Q 


(Vm.26) 


These  equatimis  combined  with 


9ik  *  --A  52^,1  +  2B{k)grk- 

I 


(VIIL27) 


^  =  JLp 

^  Mn 


A 

P  =  «(x/a)(J  +  i5;A(/)w 


(VnL28) 


yidd  in  matrix  canonical  form 

\<i] 

«•  - 

P  “ 


0  0 

0  0 

c(x/a)  b 
-Afbt  -2B 


0 

-A  +  2B 
0 
0 


(VIIL29) 


with  b  a  row  matrix  with  elements  The  matrix  A  is  the  analog  of  the  matrix  in 

(VIL29)  with  die  constant  A  arul  B  is  die  diagonal  matrix  in  (VIL30).  This  is  really  a  pardtioned 
problem,  Le, 

KI“(V  -A  +  2bI[^I 


(VIIL30) 


Ditfinendating  (VIII.30)  widi  respect  to  the  time  parameter  and  using  (V1IL31)  tme  obtains 

[<51  _f  e(,/.)/Ar.  bfM.  IfQl.pfQl  ,V1IU2^ 

This  nonhermidan  problem  can  be  brought  to  classical  cammical  form  by  a  similarity  transfor¬ 
mation  A  s  L~^FL  Le.  the  resulting  matrix  A  consists  of  Jordan  blodrs 


(VIIL31) 


(VIIL32) 


fA.  I  1  A|  1  0 

[»  A*l’  [0“  0  ij’ 


(VIIL33) 


For  die  case  where  aU  the  Jmdan  bloda  ate  tme-dimensitmal,  the  matrix  A  is  diagonal  Then, 
defining  a  new  colunm  matrix 


l8<J 


(VIIL34) 


the  equadcm  of  motion  can  be  expressed  as 

0  =  Afi.  (VIIL35) 

The  negative  of  die  eigenvalues  of  A  ate  dien  die  fiequencies  squared  of  the  coupled  system 
of  electrtms  and  nudel  The  requirement  that  the  firequoicies  be  teal  puts  furdm’  constraints 
on  the  parameters  in  the  matrices  A,  B,  and  b.  For  small  vibrations  around  the  ground  state 
an  fiequendes  should  be  real 
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DC.  Conclusions 


•  • 


The  END  the(»y,  which  provides  t  fiiUy  dynamical  tieatment  of  electrons  and  nuclei 
without  lerorting  to  the  ctmunon  Btmi-Oppenhdmer  ot  Adiabatic  approximations,  has  been 
used  to  fonnulate  an  internally  consistent  description  of  electrons  and  vibrating  nuclei  of  an 
extended  system.  In  older  to  demonstrate  the  essential  features  of  such  a  systmn,  allowing  for 
the  coupliqg  between  nuclear  and  electronic  dynamics,  a  PPP-like  hamihonian  is  mnployed. 
Such  a  hamiltnnian  yields  a  total  energy  that  scales  properly  with  the  size  of  die  system  and 
exhibits  a  jdiysicaUy  correct  dependence  on  small  nuclear  di^lacements.  It  also  satisfies  the 
{Kimary  requirement  of  all  treatments  of  extended  dectronic  systems,  namely  permitting  a  propa 
tieatment  of  the  long  range  nature  of  the  Coulomb  interaction. 

Ihe  nuclei  are  treated  classically  and  the  dectnms  are  described  by  a  single  determinantal 
wavefuncdmi  built  with  complex  spin  oibitals.  Ihe  nuclear  positions  and  conjugate  momenta  are 
die  dynamical  variables  of  the  nuclear  dynamics,  while  the  electronic  dynamical  variables  are 
the  complex  coefficients  mixing  Bloch  f^tions  below  with  those  above  the  Fdmi  level  The 
Wannier  funedons  (or  corresponding  Bloch  fiincdons)  in  one  band  associated  with  the  equilibrium 
nuctear  posidmis  "idry  up  die  electronic  basis.  The  particular  sample  ^stmn  studied  is  metallic 
with  a  half-filled  band  in  its  reference  state. 

The  equations  that  govern  the  time  evolution  of  the  dynamical  variables  are  derived  via  the 
time-dqiendent  variational  principle.  Iheir  analysis  is  presented  in  three  steps.  First,  frozen 
electronic  degrees  of  freedom  lead  to  a  variant  derivation  of  the  classical  equations  of  motion 
of  vibrating  nucld.  Second,  frozen  nuclear  degrees  of  freedom  yield  a  RPA  description  of 
the  electronic  excitations.  Ihird,  the  Adi  END  equations  are  solved  revealing  that  the  highest 
frequency  nuclear  mode  is  the  only  one  coupling  to  the  electron  dynamics. 

In  qiite  of  its  simplicity,  a  linear  chain  widi  one  atom  per  rnit  cell  and  a  single  valence  elec¬ 
tron  per  atom  exhibits  the  essential  complexities  to  demonstrate  the  generality  of  the  qiproach, 
yet  permits  a  straightforward  and  idatively  uncluttered  analysis.  Additional  specializations  in¬ 
clude  the  detailed  treatment  of  only  longitudinal  vibrations  and  the  selection  of  only  one  mixing 
coefficient  per  occupied  spin  tnbital  lii^g  it  to  a  single  virtual  ^in  orbital  as  is  done  in  the 
alternant  molecular  orbital  AMO  method. 

Permitting  a  detailed  and  transparent  application  of  the  theory  to  an  extended  system  drese 
simplificatirms  in  no  way  reflect  any  basic  limitations  of  the  END  ^iproach.  As  is  indicated 
in  sectiems  IV  and  VI  transverse  vibrations  frdlow  in  exaedy  the  same  manner  as  in  traditional 
treatments.  Similariy,  systems  with  more  than  one  atom  per  unit  cell  involve  more  work,  but  no 
additional  conceptual  challenges.  Also,  the  consideration  of  several  electronic  mixing  coefficients 
leads  to  more  dynamical  equations,  but  no  significandy  different  behavior  from  that  of  the  simple 
AMO  choice,  friclusion  of  more  than  one  energy  band  again  means  more  work  but  adds  nothing 
fundamentally  new. 

Extensions  of  the  simplest  END  model  to  one  with  quantum  nuclei  and  possibly  also  a 
multiconfiguiational  electronic  description  have  been  accomplished  for  finite  molecular  systmns 
and  might  also  be  attempted  for  extend  s^tons.  In  particular,  a  quantum  mechanical  treatment 
of  die  nuclei  is  necessary  to  permit  a  discussion  of  phoiKm  systems  and  a  study  of  electron- 
phonon  coiqiling. 
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